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Statistical mechanical models with local interactions in d > 1 dimension can be regarded as 
d = 1 dimensional models with regular long range interactions. In this paper we study the critical 
properties of Ising models having V sites, each having z randomly chosen neighbors. For z = 2 the 
model reduces to the d = 1 Ising model. For z = oo we get a mean field model. We find that for 
finite z > 2 the system has a second order phase transition characterized by a length scale L — \nV 
and mean field critical exponents that are independent of z. 



PACS numbers: 05.50, 05.70.J, 05.70.F 

The dimensionality of a system and the symmetries of 
its Hamiltonian determine its critical behavior, which is 
characterized by the exponents that describe the singu- 
larities in magnetization, susceptibility, and specific heat 
at the critical temperature T c . Meanwhile, features of 
a model such as additional terms in the Hamiltonian of 
the same dimensionality and symmetry are irrelevant to 
these properties. Such "universality" applies, however, 
only to systems with the same range of interactions. Ex- 
ponents are affected by long range forces which can alter 
the way correlations diverge at the critical temperature. 

As a particular instance of this effect, increasing the 
dimensionality can be regarded as the addition of pecu- 
liar long range interactions: The d — 2 Ising model on a 
N x N lattice is trivially equivalent to a "d — 1" model 
where site I connects to sites I ± x and I ± Ni. This new 
model is characterized by a coordination number z = 4 
and a regular set of long range bonds. The critical expo- 
nents were changed by the addition of the / ± Nx bonds. 
In order to study this effect, and in particular to under- 
stand the importance of the regularity of the additional 
connections, we consider in this paper a u d = 1" ferro- 
magnetic Ising model 

H = -JJ2S 1 S J (1) 

where Si is a classical variable which can take the values 
±1 and the sum is over "random" neighbors (ij) defined 
more precisely as follows. Suppose we have a total of V 
sites each of which is to have z randomly chosen neigh- 
bors. Our operational prescription for the construction 
of random neighbors is to begin with a list of integers of 
length zV in which the first V positive integers are re- 
peated z times. This list is randomized by a large num- 
ber of pair interchanges. Sites are defined to be neigh- 
bors if they appear as the In — 1 and In entries in the 
list, n = 1,2,...^-. We eliminate by further randomiza- 
tion any pairs that appear twice or sites that are self- 
connected. The end result of this procedure is a lattice 
in which each spin has z neighbors randomly chosen from 
the V — 1 remaining spins 0]. 



we will construct the d — 1 dimensional Ising model, if 
we insist that the lattice have all sites in the same cluster. 
For z = 4 we would have a model similar to the d = 2 
Ising Hamiltonian in the sense that each site would be 
connected to four neighbors, two of which could be ar- 
bitrarily considered as geometrically close, the other two 
of which are not. However, the disordered nature of the 
bonds has an important qualitative effect. The bonds 
are "more long ranged" than the regular d-dimcnsional 
case. Consider the average separation, {aij) c , between 
two spins, defined as the average over all pairs of spins of 
the minimum number of bonds traversed in moving from 
one to another. (The subscript "c" is included to em- 
phasize that this is an average over the different pairs of 
sites in the lattice rather than the usual thermodynamic 
average. We could also ask about an average over differ- 
ent possible realizations of the random bonds. However, 
we shall assume that in the thermodynamic limit of large 
V the lattice is self-averaging.) For short-ranged models 
in any dimension d, the average separation gr~ows with 
the linear extent L. That is, (cry) ~ L = Va. For the 
random bond case, (cry) — a(z) + b(z)hiV. The average 
separation of spins is much smaller. The role of z is not 
fundamental in the sense that it controls only the values 
of a and 6, not the logarithmic dependence on V. This 
result can be obtained numerically by explicit construc- 
tion of lattices and computing cry . In Fig. 1 we show such 
numerically obtained values for our random long range 
model. 

The behavior of the average separation with volume 
can also be motivated analytically by considering a 
Bethe lattice of coordination number z. Such a lattice 
differs from that described above in having no closed 
loops, but nevertheless, as we will see in greater de- 
tail below, is relevant to the random bond model we 
are studying. The central site from which the lattice 
is grown has z first neighbors, z(z — 1) second neigh- 
bors, and z(z — l)™ -1 neighbors of distance n. If the 
lattice is terminated after level N, then the total num- 
ber of sites is V = 1 + zJ2n=i( z ~ I)™ -1 and the dis- 
tance to the central site in the thermodynamic limit is 
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FIG. 1: The average separation between spins {a, 
as a function of \nV for different choices of the coordination 
number z. The slope of the curve is b(z). 



tice of coordination number z exhibits a similar logarith- 
mic growth of (<7y) with V and has the particular value 

This very slow growth, in a fashion qualitatively inde- 
pendent of z for z > 3, of the "distance" between bonds 
with the system size might suggest that the critical prop- 
erties will be more similar to that of the mean field case, 
where the distance between all bonds is unity, than the 
case of finite dimension where (ov,) ~ V^. Indeed, the 
apparent absence of a length scale would also suggest 
mean field like behavior. This is in accordance with the 
fact that if one defines the dimensionality of a systemQ 
as d = linin^oo [lnc n /lnn], where c„ is the number of 
sites within n steps of a given site, then one finds d = oo. 
However, the number of paths between spins is also rel- 
evant in determining the critical behavior [? ], and this 
number of paths does depend on z, so these arguments 
suggestive of mean field behavior are not conclusive. 

To make this point more precise, we note that regular 
d-dimensional lattice Hamiltonians are characterized by 
the existence of short loops of bonds which connect sites 
back to themselves. In contrast, if the underlying lattice 
is a Cayley tree, sites are connected to neighbors without 
any such closed loops. The model discussed here can be 
viewed as an alternate way of tying off the ends of a 
Cayley tree lattice in which closed loops of size InV are 
created. For such structures it is known that the manner 
in which the boundary, which is a nonvanishing fraction 
of the total number of sites in the thermodynamic limit, 
is treated is important to the critical behavior 0. 

Calculation of the energy cost for introducing a domain 
wall which divides the lattice in two is also instructive. 
Consider first the usual Ising model on a hypercubic lat- 



a single bond need be broken. In two dimensions V? 
bonds must be broken, and in d dimensions we need to 
break V~a~ bonds. This increasing difficulty of forming 
a domain wall is, of course, why the tendency towards 
ordering is higher as the dimensionality increases. In the 
random bond model we have described, an attempt to 
divide the lattice in two equal pieces would on the aver- 
age necessitate a breaking of ¥ bonds, a scaling with the 
volume V which is the d — > oo limit of the usual Ising 
case. This further points towards a mean field scenario 
for the critical properties. 

We can analyze in a little more detail the effect of 
the behavior of (<7y ) by considering a high temperature 
expansion. To lowest order 
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tanh(/?J). 



(2) 



(The Griffith's inequality in fact guarantees that this is a 
lower bound on the correlation function |5|].) In 1 dimen- 
sion for free boundary conditions the shortest path is the 
only path and this expression is exact. (SiSj) decays ex- 
ponentially to zero with increasing separation except at 
T = T c = 0. 

The magnetic susceptibility is a sum over all such cor- 
relation functions. For temperature T > T c where the 
magnetization is zero, the susceptibility per site is 
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(3) 



Thus we estimate a critical temperature at which di- 
verges in the thermodynamic limit as 



= tanh(/3 c J) e . 



(4) 



This argument, in fact, gives the exact result for the 
Bethe lattice, since if we substitute b(z) 



ln(z-l) 



into 



within MFT, 2k = z 
a plot of — 



Eq. (4), we obtain (3 C J = |ln-^-, which is precisely the 
Bethe lattice critical temperature. On the other hand, 
In Fig. 2 we show, as the full curve, 
vs. z obtained from Eq. (4) using Bethe 
lattice values of b(z) which are close to those obtained 
numerically for the random bond prescription described 
above. The squares are obtained by scaling monte carlo 
results to obtain T c as discussed shortly. 

In order to see whether these arguments suggesting 
some sort of mean field theory really are valid, we have 
performed monte carlo simulations. In Fig. 3 we show 
results for the susceptibility per site, y, vs T for z = 3 
and 4. The different symbols represent lattices of differ- 
ent size, typically between V = 200 and V = 6400 sites. 
The increasing sharpness and height with volume indi- 
cate the presence of a finite T c which increases with the 
coordination of the lattice. Fig. 4 shows similar results 
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FIG. 2: The approximate analytic form for the critical tem- 
perature T c obtained from Eq. (4) normalized by the coor- 
dination number z. The Bethe lattice form of b(z) has been 
used. The squares are values for the normalized critical tem- 
perature obtained from scaling finite size monte carlo results 
for the susceptibility and the specific heat. 

nonanalytic behavior as V — » oo. The values of T c are 
consistent with those suggested by Fig. 3. 

To explore the order of the phase transition, we can 
look at histograms of the magnetization distribution 
P{m) for different temperatures. At high temperatures 
P(m) exhibits a single peak at m — 0. As T is decreased 
this peak is broadened and eventually splits into two 
peaks symmetrically located about ±mo. At no temper- 
ature is evidence seen for the 3 peaked structure charac- 
teristic of a first order transition where m = and m/0 
phases are in coexistence. This suggests the transition is 
second order for finite z, as it must be for z = oo. 

It is natural to attempt to perform finite size scaling on 
the data. We can first imagine the maxima in the plots 
of the susceptibility and specific heat represent values of 
"T c " on finite size lattices. Conventional scaling theory^ 
suggests that 

r c (V0 = T e (oo)[l-L-"]. (5) 

We are immediately confronted with the difficulty of re- 
lating the linear extent L to the number of sites V. 
Choosing L = V would in the usual d-dimensional mod- 
els simply rescale the exponent v by the dimensionality. 
We find that to obtain a reasonable fit to this form we 
would need to use an anomalously large value of v. In 
fact, a better ansatz is obtained by replacing L by lnV. 
We then find that the best fit for z = 3 is obtained by 
T c (oo) = 1.91 ± 0.05 and v = 0.52 ± 0.06. This is, of 
course, in reasonable agreement with the mean field value 
for the correlation length exponent v. 
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FIG. 3: Raw results for the susceptibility per site as a func- 
tion of temperature for different sized lattices. Fig. 3a is for 
z = 3 and Fig. 3b for z = 4. 

X, we expect that curves for different lattices, distinct in 
their unsealed forms of Figs. 3 and 4, will fall on univer- 
sal curves when scaled with appropriate choices of the 
critical temperature and exponents 6]. That is, for \ 

jj =f(Lh), (6) 

while for C 

^ ¥ =g(Lh). (7) 

Here / and g are some scaling functions. In Fig. 5 we 
show x scaled appropriately using L=lnV* and the crit- 
ical exponents v = 0.53 and 7 = 1.06. We have chosen 
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FIG. 4: Raw results for the specific heat per site as a function 
of temperature for different sized lattices. Fig. 4a is for z — 3 
and Fig. 4b for 2 = 4. 



The critical temperature T c = 1.87. This scaling plot 
has somewhat more scatter than the ones obtained for 
regular models 0- However, this can be attributed to 
the need to average the results, especially for smaller lat- 
tices, over different configurations of random bonds. For 
the specific heat, we note that the maximum in C in 
Fig. 4 is insensitive to lattice size. This translates to a 
small value for the exponent a, which is zero in the mean 
field limit. To estimate the accuracy of these values of 
the exponents, we have examined scaling plots in which 
the exponents have been shifted. We see a significant 
deterioration of the quality for changes of more than 10 
percent. 

We have seen that the Ising model defined by Eq. 1 and 



o 



X 



ro - 



CM 



O 



T c = 1.870 
z=3 



4- ± 



-15 -10 -5 5 
(T-T c ) ln(V) b 

FIG. 5: The results of Fig. 3a are shown using the appro- 
priate scaled forms as described in the text. The critical ex- 
ponents used were v = 0.53 and 7 = 1.06. The critical tem- 
perature T c =1.87. Different symbols correspond to different 
lattice sizes, K=800, 1600, 3200, and 6400. 



field like critical behavior with a length scale which de- 
pends logarithmically on the "volume" of the system. 
Clearly the regular d-dimensional models are realized in 
our model by certain (very) improbable choices of the 
random bonds. We can ask about the renormalization 
group flows in the space of possible choices of the bonds. 
If we construct some quantitative measure of the distance 
of a random assignment from the usual regular one, what 
will be the extent of the basin of attraction of the fixed 
points corresponding to the regular models and their ex- 
ponents? Is it vanishingly small, so that in some sense all 
models of this sort flow towards mean field theory, or will 
there be some finite region characterized by the regular 
exponents? 

Traditionally, one route to mean field theory is to in- 
crease the dimensionality of the system, which has the 
effect of enhancing the connectivity of the lattice while 
leaving the interactions short ranged. In a similar man- 
ner, mean field results are obtained for the critical prop- 
erties of the Bethe lattice in the limit where the coor- 
dination z goes to infinity, but again the interaction is 
short ranged. Both these approaches, then, take a route 
to MFT via increased numbers of neighbors. However, 
we have seen here a novel realization which does not pro- 
ceed via increased connectivity but rather works with 
fixed coordination and considers a model with extremely 
long ranged interactions. Kac et.al.yj have shown in a 
one dimensional model with an interaction of length scale 
R that there is a nonzero T c only in the R — » 00 limit. 
There, however, R and z are simultaneously increased. It 
is clear from our work that increasing R alone, at fixed z, 
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that there are two alternate routes to mean field critical 
behavior: cither diverging z or diverging R alone suffice. 
As we have pointed out, this is perhaps not too surpris- 
ing since neighbors in additional spatial dimensions can 



equivalently be regarded as regular long range bonds. 
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